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I.  INTRODUCTION  AND  SUMMARY  OF  THE  RESEARCH  PROGRAM 


To  properly  describe  coupled 
chemical  reactions  and  gaseous 
diffusion  in  porous  sorbent  and 
catalyst  grains,  the  "pore  tree"  was 
introduced  by  Simons  and  Finson 
(1979)  and  Simons  (1982).  The  pore 
tree  represents  an  isolated  sub¬ 
structure,  allowing  diffusion  into  and 
out  of  porous  media  without 
permitting  transport  through  the 
media.  A  pore  size  distribution  was 
derived  statistically  and  confirmed 
empirically  for  coal,  coal  char, 
sorbents,  catalysts  and  kidney  stones 
from  both  men  and  women.  The 
pore  tree  was  derived  from  the  pore 
size  distribution  and  allows  the 


Figure  1.  The  Pore  Tree 


orderly  migration  of  a  reactant  gas 

from  die  large  pores  to  the  smaU  pores  (Fig.  1).  A  detailed  description  of  the  pore  tree  and 
the  coupled  transport  and  chemistry  is  given  by  Simons  (1982,  1983a).  The  spatially 
dependent  transport/reaction  equations  are  solved  for  a  single  pore  tree  and  then  the  total 
contribution  of  all  trees  (of  all  sizes)  in  the  system  is  obtained  by  summing  the  contribution  of 
each  tree  that  reaches  the  exterior  of  the  system.  This  is  distinct  from  the  "bulk"  transport 
approach  in  which  the  transport  equation  for  a  single  pore  is  integrated  over  all  pores  at  a 
fixed  point  in  space  before  integrating  spatially.  The  "bulk"  transport  approach  is  invalid  if 
the  spatial  gradients  in  the  transport  equations  are  implicit  functions  of  pore  size.  One 
example  of  this  implicit  pore  size  dependence  is  that  of  the  heterogeneous  reactions  within 
porous  sorbents  (Simons,  1988)  and  catalysts  for  which  the  pore  tree  structure/transport 
model  was  developed.  A  second  example  is  that  of  coupled  diffusion  and  remediation 
reactions  in  the  immobile  region  of  soil. 


In  order  to  describe  the  subsurface  transport  of  gas  and  water  in  soU,  the  dispersion  of 
contanunants,  and  in-situ  remediation  of  contaminated  sites,  the  pore  tree  has  been  extended 
(Simons,  1996a)  to  simulate  permeability  and  bulk  transport.  The  interconnectivity  of  the  pore 
structure  is  illustrated  in  Fig.  2  and  is  obtained  via  a  statistical  determination  of  the 
"branches"  that  are  common  to  several  trees  to  allow  convection  and  bulk  diffusion  through 
the  large  scale  (mobile)  structure  in  addition  to  diffusion  and  coupled  chemical  reactions 
within  the  smaller  scale  (immobile)  structure.  The  statistical  analysis  has  determined  that  the 
probability  of  pore  interconnectivity  extends  across  the  entire  pore  size  range,  with  an 
increase  in  the  probability  accompanying  a  decreasing  pore  size.  While  permeability  is 
dominated  by  the  largest  pores,  it  is  also  important  to  establish  the  level  of  convection  and 
diffusion  that  is  occurring  at  the  intermediate  scales  in  order  to  accurately  relate  large  scale 
bulk  transport,  small  scale  diffusion  and  coupled  chemical  reactions. 
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a)  Pores  Interconnected  in  Plane  AA  b)  Pores  Interconnected  Out  of  Plane  AA 

Figure  2.  Interconnectivity  of  the  Pore  Tree 


The  permeability  across  a  given 
plane  is  limited  by  the  largest  pores  that  are 
interconnected  in  that  plane  (Fig.  2a).  The 
statistical  analysis  has  determined  that 
approximately  one  quarter  of  one  percent  of 
all  large  pores  are  interconnected.  This 
establishes  a  very  coarse  grid  for  the 
permeability  which  leads  to  the 
measurement  scale  size  errors  illustrated  in 
Fig.  3.  The  extended  pore  tree  model  has 
successfully  explained  the  measurement 
errors  in  the  permeability  of  soil  due  to  the 
measurement  scale  size  (Shouse,  et.al., 

1994)  which  has  indirectly  confirmed  the 
low  probability  of  the  interconnectivity. 

The  btilk  gaseous  diffusivity  across  a 
given  plane  is  shown  to  be  limited  by  the 
interconnectivity  of  the  smaller  branches 
outside  of  that  plane  (Fig.  2b).  These  small 
pores  may  be  saturated  and  gas  diffusion  is 
allowed  only  in  pores  whose  radius  is 
greater  than  that  of  the  largest  saturated 
pore  (r,g^.  The  gas  diffusivity  is  shown  to 

scale  as  (permeability) r^.  A  comparison  of  the  present  theory  to  the  diffusivity  data  of 
Washington  eL  al.,  (1994)  is  illustrated  in  Fig.  4.  The  gas  diffusivity  data  at  approximately 
60%  saturation  verifies  the  permeability  dependence  and  suggests  a  saturation  r^us  of  the 
order  of  30  pm.  Sixty  percent  saturation  at  30  pm  is  consistent  with  the  statistically  derived 


Figure  3.  Measurement  Scale  Errors 
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CONDUCTIVITY  (cnVhr) 


Figure  4.  Bulk  Diffiisivity  of  a  Partially  Saturated  Soil 

pore  size  distribution  function.  While  the  excellent  agreement  with  the  data  does  substantiate 
the  present  theory,  the  diffiisivity  in  partially  saturated  soU  is  very  sensitive  to  an  unknown 
saturation  radius.  If  bulk  diffusivities  are  to  be  correlated  with  field  data,  such  measurements 
should  attempt  to  measure  the  saturation  radius. 

The  permeability  and  the  bulk  diffiisivity  have  tested  two  extreme  limits  of  the  pore 
structure  and  pore  interconnectivity  concepts.  Permeability  is  limited  by  the  in  plane 
interconnectivity  (Fig.  2a)  and  bulk  gaseous  diffusion  is  limited  by  the  out  of  plane  (Fig.  2b) 
interconnectivity.  Permeability  is  limited  by  the  large  pore  interconnectivity  and  bulk  diffusion 
is  limited  by  the  interconnectivity  of  the  smaller  pores.  The  apparent  success  of  these 
concepts  over  a  very  broad  pore  size  range  suggests  that  the  extended  pore  tree  model  will 
accurately  describe  the  relationship  between  large  scale  convection  and  small  scale  diffusive 
transport. 

A  preliminary  step  in  this  approach  to  relate  convective  and  diffusive  transport  in 
various  size  pores  is  to  describe  the  subscale  convection  responsible  for  hydrodynamic 
dispersion  (Simons,  1996b).  The  velocity  profile  in  interconnected  pore  space  is  illustrated  in 
Fig.  5.  The  largest  pores  are  not  interconneced  across  the  porous  medium.  Hence,  the  net 
velocity  within  the  largest  pores  is  zero.  The  smallest  pores  do  not  support  a  very  large 
convective  velocity.  Thus,  there  is  a  relative  maximum  of  the  velocity  at  some  intermediate 

pore  size  .  Clearly,  the  fluid  in  the  pore  of  radius  convects  ahead  of  that 

in  the  pores  of  radius  ty  and  X2.  Consider  flow  in  a  saturated  porous  medium  in  which  all  fluid 
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Figure  5.  Velocity  in  Pore  Space 
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Figure  7.  Temporal  Concentration  Profile 
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Figure  6.  Spatial  Concentration  Profile 


Figure  8.  Temporal  Concentration  Data 


behind  a  fixed  plane  normal  to  the  flow  is  suddenly  injected  with  red  dye.  If  fluid  velocity 
v„*  is  associated  with  pores  ri*  and  rj*  ,  aU  dyed  fluid  in  pores  whose  radius  is  between  rj* 
and  r2*  will  translate  with  velocity  equal  to  or  greater  than  v„*  and  move  to  a  corresponding 
location  X*  or  greater.  The  corresponding  concentration  profile  is  illustrated  in  Fig.  6  where 
all  red  dye  would  have  reached  the  location  Eta  =  1  were  it  not  for  the  velocity  variation 
within  pore  size.  The  "front"  of  the  concentration  profile  has  progressed  30%  ahead  of  the 
mean  while  the  slower  fluid  in  the  smaller  pores  has  trailed  considerable.  If  we  were  standing 
at  the  end  of  a  porous  flow  tube,  the  concentration  passing  our  location  would  exhibit  a  time 
dependence  illustrated  in  Fig.  7.  The  faster  fluid  breaks  through  ahead  of  the  mean  flow 
while  the  fluid  in  the  smaller  pores  is  delayed. 

Data  of  EWch,  et  al.,  (1966),  as  reported  by  Brusseau  and  Rao  (1989),  are  illustrated 
in  Fig.  8  and  support  the  current  theory  of  hydrodynamic  dispersion.  While  the  long  "tail"  of 
the  concentration  profile  is  generally  attributed  to  nonequUibrium,  the  current  pore  structure 
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model  attributes  it  to  the  wide  pore  size  distribution  occurring  in  soil.  The  narrow  pore  size 
distributions  occurring  in  laboratory  "soil”  will  yield  a  much  narrower  concentration  profile, 
i.e.,  less  hydrodynamic  dispersion.  Future  research  will  be  directed  toward  1),  evaluating  the 
hydrodynamic  dispersion  as  a  function  of  pore  size  range  in  the  soil  sample  and  2),  using  this 
analytic  pore  structure/  pore  transport  model  to  help  derive  an  analytic  expression  for 
hydrodynamic  dispersion  that  replaces  the  Fickian  process  in  the  transport  equations. 

Analysis  of  the  permeability,  bulk  diffusivity,  and  hydrodynamic  dispersion  has 
utilized  the  interconnectivity  of  the  pores  to  determine  the  distribution  of  the  convection 
velocity  with  pore  size.  This  analysis  suggests  the  existence  of  a  permeable  sub-range  in  the 
pore  structure  which  does  not  contribute  significantly  to  the  bulk  permeability  but  in  which 
convection  dominates  diffusion.  It  is  the  bidance  of  the  sub-scale  convection  with  the  small 
scale  diffusion  that  will  control  contaminant  transport  and  in-situ  remediation.  The  size 
distribution  of  the  pores  and  grains,  and  the  variations  in  fluid  velocity  within  and  between 
pores  of  different  sizes  is  critical  to  interfacing  the  transport  processes.  A  methodology  has 
been  developed  to  couple  subscale  diffusion,  convection  and  chemical  reactions  to  the 
macroscopic  transport  in  order  to  accurately  describe  contaminant  transport  and  in-situ 
remediation  in  Ground  Water  Simulation  codes. 

A  typical  species  transport  equation  in  permeable  soil  is  of  the  form: 

dc  -dc  ^  ^  cM 

dt  ~M 

where  c  is  the  local  species  concentration,  v  is  the  mean  convection  velocity  as  determined 
by  the  local  pressure  gradient  and  permeability  k ,  and  is  the  hydrodynamic  dispersion. 
The  dcldt  term  is  the  true  unsteady  term  and  the  source  and/or  sink  of  species  c  due  to 
chemical  and/or  physical  processes  within  the  subscale  pore  structure  is  written  as  cMjM 
where  M/M  represents  a  bulk  rate  (1/t)  of  production  or  consumption.  The  bulk  rate  M/M 
will  be  related  to  gradients  of  c  on  length  scales  of  order  millimeters  within  the  subscale  pore 
structure  and  cannot  be  expressed  in  terms  of  the  dcjdx  of  the  macroscopic  transport  grid. 

To  describe  the  subscale  chemical  reactions  and  transport,  a  "grain"  of  soil  of  radius 
Ug  is  isolated  from  the  rest  of  the  medium.  The  grain  size  will  be  chosen  sufficiently  small 
that  diffusion  and  coupled  chemical  reactions  will  dominate  the  subscale  transport.  The  pore 
tree  was  developed  (Simons,  1982,  1983a,)  to  treat  coupled  chemical  reactions  and  diffusion 
within  immobile  porous  grains.  Generic  forms  of  the  immobile  solution  are  illustrated  in 

Figure  9  where  /?,  represents  the  limit  of  kinetic  control  and  represents  a  coupled 
kinetic/diffusive  solution  in  which  the  reaction  on  the  walls  of  the  pore  establishes  a  species 
gradient  which  accelerates  the  molecular  diffusion.  The  particular  solution  for  contaminant 
transport  and  soil  remediation  will  vary  with  the  number  of  reactions  and  kinetic  mechanisms 
but  the  generic  form  of  the  immobile  solution  will  be  similar  to  that  illustrated  in  Fig.  9. 
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Figure  9.  Coupling  The  Mobile  and  Immobile  Solutions 


The  reactant  species  must  be  supplied  to  the  edge  of  the  grain  by  subscale  convection 
at  a  rate  depicted  by  the  permeability  which  is  appropriate  to  the  length  scale  ag.  The  mobile 

solution  for  M/Af  is  developed  in  subsequent  sections.  It  is  shown  that 

^  ^  <^ldx)  s  {Convective Limit) 

^  768ppA:^  * 


Tiie  mobile  solution  indicates  that  the  convective  limit  increases  with  grain  size  while 
the  immobile  solution  indicates  that  the  diffusive  limit  decreases  with  increasing  grain  size. 
The  comparison  illustrated  in  Figme  9  suggests  that  for  large  grains,  diffusion  cannot  keep  up 
with  convection  while  for  small  grains,  convection  cannot  keep  up  with  diffusion.  The 
convective  and  diffusive  limits  are  in  balance  if  and  only  if 

^-SmaawQf 


where  it  has  been  noted  that  the  solution  cannot  exceed  the  limit  of  kinetic  control. 

The  above  expression  for  the  bulk  source  and/or  sink  of  species  c  due  to  chemical 
and/or  physical  processes  within  the  subscale  pore  structure  is  incomplete  without  specific 

models  for  and  in  the  immobile  region.  Future  research  will  be  directed  toward 
developing  a  library  of  such  models.  The  most  basic  processes  of  contaminant  diffusion  and 
adsorption  are  briefly  described  below. 
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n.  ISOLATED  PORE  TREE:  THE  STRUCTURE 


Following  the  pore  structure  theory  of  Simons  and  Finson  (1979)  and  Simons  (1982), 
consider  a  spherical  porous  particle  of  radius  a,  containing  pores  of  length  Ip  and  radius  rp. 
The  pore  dimensions  range  from  a  microscale  of  the  order  of  Angstroms  to  a  macroscale 
which  is  a  significant  fraction  of  the  particle  radius.  The  radius  of  the  largest  pore  is  denoted 
by  r^  and  is  given  by 


'•„=2a0‘P/3K'^  (1) 

where  6  is  the  total  porosity  of  the  particle  and  is  a  constant  of  integration,  approximately 
equal  to  five,  which  relates  the  pore  length  to  its  radius 

The  radius  of  the  smallest  pore  is  denoted  by  r^  and  is  given  by 


where  p,  is  the  density  of  the  solid  matrix,  Sp  is  the  specific  internal  surface  area  (several 
hundred  m7g),  and 


P  =Mr^/r^) 


(4) 


The  particle  contains  a  continuous  distribution  of  pore  sizes  fix)m  r^^i^  to  r,„g„  .  The 
number  of  pores  within  an  arbitrary  plane  of  cross-sectional  area  A  and  with  radius  between 
rp  and  rp  -hdrp  is  denoted  by  girp)Adr^.  The  pore  distribution  function  5(r^)is  given  by 


g(rp=0/2npr;  (5) 

where  indicates  an  average  over  all  inclination  angles  between  the  axis  of  the  pore  and 
the  normal  to  the  plane.  Due  to  the  random  orientation  of  the  pores,  the  intersection  of  a 
circular  cylinder  with  a  plane  is  an  ellipse  of  average  area  27crp^.  Hence,  the  porosity  is  the 
2jirp^  moment  of  g(r^)  and  the  internal  surface  area  is  the  47irp  moment  of  g(r^). 

The  above  expression  for  was  derived  firom  statistical  arguments  and  was 
subsequently  validated  for  coal  char  through  comparison  of  the  predicted  volume  and  surface 
area  distributions  with  mercury  intrusion  data  (Simons  and  Finson,  1979).  Mercury  intrusion 
data  generally  demonstrate  a  linear  increase  in  the  intrusion  volume  with  In(rp).  It  is  the 
functional  form  of  this  relationship. 
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(6) 


y 

Pore  Volume  «  j rpg{r^)dr^ 


'min 


that  depicts  the  inverse  cubic  dependence  of  g(r^)  on  fp.  Mercury  intrusion  data  has  been 
used  to  validate  this  pore  size  distribution  function  for  coal  and  char  derived  from  that  coal 
(Kothandaraman  et.al.,  1984),  sorbents  (Simons  and  Garman,  1986),  catalysts,  and  even 
kidney  stones  from  both  men  and  women.  There  is,  however,  no  such  validation  for  soil  or 
sediments  available. 

The  number  of  pores  within  the  bulk  volume  V  whose  pore  radius  is  between  rp  and 
rp+drp  may  be  defined  by  K/(r^)  dr^.  The  pore  volume  is  expressed  as  the  Jcrp^lp  moment  of 
/(r^)  and  the  internal  surface  area  is  the  27irplp  moment  of  /(r^) .  The  pore  size  distribution 
functions  ( f(r^)  and  g(r^)  )  are  clearly  not  independent.  The  definitions  of  porosity  and 
internal  surface  area  infer  that  f(r^)  is  related  to  gir^)  by 


gifp)  =f(rp)lpl2 


(7) 


Equation  (7)  simply  states  that  the  probable  number  of  pores  intersecting  an  arbitrary  plane 
increases  with  the  length  of  the  pore  and  with  the  density  of  pores. 


The  length  of  a  pore  is  determined  by  an  arbitrary  intersection  with  another  pore  and 
is  expressed  (Simons  and  Finson,  1979)  as  a  collision  integral  over  the  pore  distribution 

functions.  The  analysis  suggests  that  Ip  ,  J(r^)  and  /(r^)  are  proportional  to  rp ,  Ifrp^  and 
l/rp‘*  respectively.  The  constants  of  proportionality  are  obtained  from  integral  constraints,  i.e., 
the  total  porosity  and  internal  surface  area  contained  in  the  pore  structure.  I'he  expression  for 
fitp)  is  given  by 


f(rp)  = 


Q*P 


(8) 


where  the  constants  were  defined  above. 

The  pore  volume  distribution  corresponding  to  these  distribution  functions  is  similar  to 
that  utilized  in  the  random  pore  model  (Gavalas,  1980  &  1981).  However,  the  pore  tree 
model  and  the  random  pore  model  differ  dramatically  in  their  choice  of  the  pore  aspect  ratio 
(length  to  diameter)  and  its  implications  with  respect  to  pore  branching.  The  random  pore 
model  allows  a  single  pore  to  connect  two  larger  pores.  This  picture  lends  itself  to  the 
idealization  of  instantaneous  mixing  between  the  pores  and  requires  that  the  pore  aspect  ratio 
be  of  the  order  of  one  hundred.  The  pore  tree  theory  uses  data  for  to  imply  (via  K^)  that 
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all  pores  possess  an  aspect  ratio  of  the  order  of  ten.  Hence,  small  pores  may  connect  to  larger 
pores  only  on  one  end  and  all  pores  must  branch  from  successively  larger  pores  like  a  tree  or 
river  system. 

Each  pore  that  reaches  the  exterior  surface  of  the  particle  is  depicted  as  the  trunk  of  a 
tree.  The  size  distribution  of  tree  trunks  on  the  exterior  surface  of  the  particle  is  denoted  by 

g(r^)4na^dr^  where  g(r^)  is  functionally  identical  to  Each  trunk  of  radius  r^  is 

associated  with  a  specific  tree-like  structure.  Let  N,  be  defined  as  the  branch  distribution 
function  where  N,drp  is  the  number  of  pores  of  radius  rp  (within  size  range  dTp)  in  a  tree 
whose  trunk  radius  is  .  The  total  number  of  pores  of  radius  rp  in  a  sphere  of  radius  a 

may  be  expressed  as  4l3na^fir^)dr^  or,  as  the  sum  of  all  pores  of  radius  rp  contained 
within  every  tree  in  the  porous  sample,  plus  all  pores  of  radius  rp  that  are  themselves  the 
trunk  of  a  tree.  Hence, 

=  f  Nfgir,)4%a^dr^  +  4na^g(r^)  (9) 


where  only  trees  with  trunk  radius  greater  than  rp  will  contain  a  pore  of  radius  rp .  Using  the 
previously  derived  expressions  for  r^,  ^(r^)and  /(r^),  Eq.(9)  is  identically  satisfied  by 


(10) 


The  branch  distribution  function  completely  characterizes  the  pore  tree.  The  internal 
surface  area  and  pore  volume  associated  with  each  pore  tree  are  denoted  by  St(rd  and  V,(rd, 
respectively,  and  are  expressed  as  the  sum  of  the  contributions  from  the  trunk  and  that  from 
the  branches. 


•t 

^nbi 

Using  Eq.(lO)  for  N„  St(rJ  and  Vt(r^  become 


(11) 


(12) 
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(13) 


S,(r,)  =27ir,/, 


<  r.  ^ 


'min 


(1-0) 


(14) 


where  the  (1-6)  term  in  St  has  been  included  to  account  for  pore  combination  (Simons, 
1979a). 


The  surface  area  associated  with  the  pore  tree  may  be  several  orders  of  magnitude 
greater  than  the  surface  area  of  the  trunk.  However,  the  volume  of  the  pore  tree  may,  at 
most,  be  one  order  of  magnitude  greater  than  that  of  the  trunk.  It  should  also  be  noted  that 
the  above  expressions  for  S,  and  reduce  to  those  appropriate  to  a  single  cylindrical  jwre  in 
the  limit  of  (the  leaf  of  the  tree).  Furthermore,  the  integrals  of  S^rJ  and  V,(r^  over 

all  g(rj) 'recover  the  total  internal  surface  area  and  pore  volume  of  the  porous  sample. 

Each  trunk  of  radius  rj  is  associated  with  a  specific  tree-like  structure  with  continuous 
branching  to  ever  decreasing  pore  radii.  The  radius  and  number  of  pores  is  a  unique  function 
of  the  distance  x  into  the  tree.  The  coordinate  x  is  skewed  in  that  it  follows  a  tortuous  path 
through  the  branches  of  the  tree.  Let  n(x)  represent  the  number  of  pores  of  radius  at 
location  X  in  a  tree  of  trunk  radius  q.  An  analysis  (Simons,  1982)  of  this  pore  tree  has 
demonstrated  that 


nix)=rflr^ix)  (15) 

and  the  coordinate  x  is  related  to  rp  by 

(16) 


The  continuous  branching  model  has  been  used  to  successfully  describe  char  oxidation 
(Simons,  1979b;  Lewis  and  Simons,  1979;  Simons,  1982  &  1983a),  coal  pyrolysis  (Simons, 
1983b  &  1984)  and  the  catalytic  cracking  of  benzene  by  porous  iron  oxides  (Simons  et  al., 
1986).  It  was  also  used  to  successfully  describe  sulfur  sorption  (SOj  and  H2S)  by  porous 
calcine  (CaO)  in  the  limit  of  zero  utilization  (Simons  and  Rawlins,  1980;  Simons  et  al,,  1984) 
and  was  later  extended  to  include  CaS04  and  CaS  deposits  (Simons  and  Garman,  1986; 
Simons  et  al.,  1987;  Simons,  1988).  The  subsequent  determination  of  the  controlling  physical 
parameters  led  to  a  new  concept  for  the  optimization  of  the  sulfur  sorption  process  (Simons, 
1991;  Simons  et  al.,  1992)  through  spray  drying  of  water  soluble  organic  calcium  solutions  to 
control  the  sorbent  pore  structure. 
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m.  INTERCONNECTTVITY 


The  first  step  in  determining  the  size  distribution  of  the  interconnected  pores  and  the 
distribution  of  the  permeability  is  to  determine  the  distribution  function  Gf(r^,r^dr^  which 
represents  the  number  of  pores  of  radius  rp  (within  size  range  drp)  per  unit  cross  section  of  an 
arbitrary  plane  and  also  contained  within  a  tree  whose  trunk  radius  is  r,.  Consider  an  infinite 
homogeneous  isotropic  porous  medium  and  isolate  a  spherical  volume  of  that  medium  denoted 
by  the  radius  a.  Such  a  volume  is  illustrated  in  Fig.  10.  The  total  number  of  pores  of  radius 
rp  (within  size  range  drp)  intersecting  plane  AA  of  area  Tca^  has  previously  been  defined  by 

g{r^)%a?'dr^  .  The  pores  in  plane  AA  in  this  size  range  may  also  be  determined  by 

integrating  G,(r^,rpiia^dr^  over  all  trees  whose  trunk  intersects  the  exterior  surface  of  the 
porous  sample.  Hence  it  follows  that 


g{rp)‘Ra^dTp  =  J  [G,(r„rpiia^dr^]  (17) 

where  only  those  trees  whose  trunk  radius  is  greater  than  rp  may  contain  a  pore  of  radius  rp . 

A  solution  to  Eq.  (17)  for  G,(r^,rp  will  not  necessarily  be  unique.  Physical  arguments 

will  help  determine  G,(r,,rp  and  help  ensure  that  it  is  the  particular  solution  we  seek.  Since 
N,  represents  the  number  of  pores  of  size  rp  in  the  tree  and  the  probability  of  a  pore 

intersecting  a  plane  is  proportional  to  its  length,  it  follows  that  should  be 


Figure  10.  Spherical  Volume  of  a  Porous  Medium 
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proportional  to  the  product  of  N,  and  lp/1, ,  i.e,,  proportional  to  r,‘/rp^,  Eq.  (17)  is  identically 
satisfied  by  a  function  which  differs  from  r,"/rp^  by  In(rp). 


4iia^r-ln(r 


(18) 


Note  that  ln(rp/r^)  introduces  an  integrable  singularity  at  rp=r„^  such  that  G,ir^,r^)  dr^  is 
finite  at  rp=rpujx  .  Hence,  there  is  one  and  only  one  largest  pore  for  each  reference  sphere. 

The  probability  of  trees  sharing  common  branches,  i.e.,  the  interconnectivity  of  the 

pore  structure  is  described  in  Fig.  11.  We  seek  the  distribution  function  which 

represents  the  number  of  pores  (within  size  range  drp  about  tp)  per  unit  area  of  plane  AA  that 
are  connected  to  both  sides  of  the  pore  structure  through  pores  at  least  as  large  as  rp.  A^  is 
defined  as  the  area  within  plane  AA  that  is  opea  to  one  side  of  the  porous  medium  through 
all  trees  of  size  r,'  (through  all  pores  of  size  rp’  that  are  at  least  as  large  as  rp).  Subsequently, 

represents  the  number  of  pores  of  size  rp  (within  size  range  drp)  per  unit  area 
of  plane  AA  that  are  contained  in  a  tree  of  size  range  dq  about  r,  and  are  also  connected  to 
the  opposite  side  of  the  porous  medium  through  all  trees  denoted  by  r,'.  It  follows  that  the 
distribution  function  for  interconnected  pores  in  plane  AA  may  be  obtained  by  integrating 

over  all  trees  (rj  that  are  large  enough  to  contain  a  pore  of  size  rp .  Hence, 


/(r^Tifl^dTp  =  J  [Ag  J(r,)27iflVr,  (19) 


Figure  1 1 .  Interconnectivity  of  a  Porous  Medium 
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From  the  above  definition  of  ,  Ag  may  be  expressed  as 


'mw  'r 

f  if  rp </rp  g(r^) 2 n a^dr^  (20) 


where  the  primes  on  the  variables  of  integration  have  been  omitted.  Evaluating  Eq.  (20) 
yields 

_iifl^01n(r^/rp 

e  2p 

from  which  Eq.  (19)  yields  the  common  branch  distribution  function. 


4P 


(22) 


It  has  been  deduced  that  the  total  number  of  common  branches  of  size  rp  in  an 
arbitrary  plane  scales  approximately  with  the  total  number  of  pores  of  that  size  in  that  plane. 
Hence,  there  is  a  probability  of  interconnectivity  at  pore  size  rp  that  is  logarithmic  in  pore 
size.  Defining  this  probability  as  Pi(rp)  via  Eq.  (22), 


= 


4p 


(23) 


it  is  apparent  that  approximately  one  percent  of  all  pores  of  all  sizes  are  interconnected 
through  larger  pores. 

The  broad  size  range  associated  with  the  interconnectivity  suggests  that  a  very  wide 
range  of  pore  sizes  control  transport  and  that  a  complicated  mixture  of  convective  and 
diffusive  transport  persists  through  all  of  pore  space.  While  permeability  is  dominated  by  the 
largest  pores,  it  is  important  to  determine  the  level  of  convection  that  is  occurring  in  smaller 
pores  in  order  to  accurately  describe  the  fine  scale  transport  necessary  to  assess  chemical 
reactions.  The  ability  of  this  pore  structure  model  to  describe  bulk  permeability,  bulk 
diffusivity  and  the  small  scale  convection  responsible  for  hydrodynamic  dispersion 
(Garabedian  et  al.,  1991)  are  demonstrated  before  attempting  to  couple  the  chemical  reactions 
to  the  macroscopic  transport. 
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IV.  PERMEABILITY 

Poiseuille's  law  for  laminar  flow  through  a  cylindrical  pore  relates  the  volume  flow 
rate  Q  to  the  pore  radius  rp,  the  viscosity  of  the  fluid  p,  and  the  pressure  gradient  — 


(f 

^  8  p  \  dtc 


whereas  the  average  volume  flow  rate  Q  across  the  cross  sectional-  area  A  of  a  porous 
medium  is  used  to  define  the  bulk  permeability  (k)  by  Darcy's  law 


A{-dpldx) 


To  determine  the  permeability,  the  pore  volume  flow  rate  must  be  related  to  the 
average  volume  flow  rate.  Convection  across  plane  AA  in  Fig.  2  will  possess  contributions 
from  two  primary  sources  illustrated  in  Figs.  2a  and  2b.  Fig.  2a  illustrates  the  case  where  the 
convection  in  plane  AA  is  due  solely  to  the  pores  that  are  interconnected  in  that  plane. 

Fig.  2b  illustrates  the  case  where  the  convection  in  plane  AA  is  due  to  the  smaller  pores  in 
the  pore  tree  that  are  interconnected  outside  of  plane  AA.  This  connectivity  will  translate  into 
a  slower  velocity  in  the  pore  crossing  plane  AA  but  could  be  significant  because  99%  of  the 
pores  in  plane  AA  are  not  interconnected  in  that  plane. 

Consider  any  pore  of  radius  r,  in  plane  AA  of  Fig.  2b  to  be  the  trunk  of  a  tree.  Each 
pore  of  size  rp  within  the  tree  possesses  the  probability  Pjir^)  of  being  interconnected  and 

each  interconnected  pore  in  the  tree  will  carry  volume  flow  rate  Qp(r^).  Since  there  are 
N,drp  (Eq.  10;  N,=r,^/rp‘*)  pores  in  size  range  drp  within  the  tree,  the  total  volume  flow  rate 

Q^ir)  through  trunk  r,  in  plane  A  A  becomes 


i>,(r,)  N,  dr^ 


or,  for  r^>  >  r.  , 


QM 


ettr^lnCr^/rJ/ 


32pp 


Higher  Order  Terms 
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Within  this  approximation,  it  is  seen  that  is  identical  to  the  volume  flowing 

through  the  pores  that  are  interconnected  within  plane  AA.  i.e., 

<?.(<•,)  ^  <?,(>•,)  W  <28) 

which  demonstrates  that  all  volume  flow  through  plane  AA  in  pore  size  r,  is  dominated  by  the 
interconnectivity  of  size  r,  in  plane  AA  and  not  by  the  interconnectivity  of  smaller  pores  in 
subsequent  branches  of  the  pore  tree.  Simply  stated:  case  2a  dominates  case  2b. 

Since  all  volume  flow  through  plane  AA  is  limited  by  the  interconnectivity  of  the 

pores  in  that  plane,  the  total  volume  flow  rate  Q  is  obtained  by  integrating  Eq.  (24)  over  all 
interconnected  pores  in  area  A.  Hence, 


(29) 


where  /(r^)  is  the  "common  branch  distribution  function"  given  by  Eq.  (22).  The  bulk 
permeability  (k)  is  then  expressed  as 


k  = 


(30) 


Upon  integration,  Eq.  (30)  becomes 


r  (31) 


Equation  (31)  resembles  a  dozen  other  expressions  (Dullien,  1979)  for  permeability 
wherein  it  is  concurred  that  the  bulk  permeability  is  dominated  by  the  largest  pores  in  the 
medium  but  the  unknown  value  of  that  permeab^ty  is  simply  replaced  by  an  unknown  fiber 
or  pore  size  to  the  second  power.  Since  the  pore  size  distribution  function  will  be  least 
accurate  at  the  extreme  end  of  the  size  range,  i.e.  at  r^  ,  no  claim  can  possibly  be  made  that 
the  numerical  constants  in  Eq.  (31)  are  in  any  way  superior  to  those  derived  elsewhere.  One 
important  advantage  of  the  extended  pore  tree  model  is  that  it  characterizes  the  distribution  of 
permeability  in  pore  space,  a  feature  that  will  be  important  in  describing  fine  scale 
contaminant  transport  and  in-situ  remediation.  A  second  advantage  is  the  ability  to  assess 
statistical  errors  in  the  measurement  of  the  permeability  as  a  function  of  the  measurement 
scale  size.  This  exercise  is  also  a  good  test  of  the  extended  pore  tree  model. 
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Consider  a  soil  sample  with 

Conductivity: 

Permeability: 

Porosity: 

Pore  Aspect  Ratio: 
ln(r^/r„J 

From  Eq.  (31),  it  follows  that 
r^  =  300  pm 

and  subsequently  r,„i„  =  20  A.  The  size  of  the  smallest  pore  is  not  an  important  parameter  for 
this  application  but  may  be  readily  adjusted  through  a  minor  variation  in  the  value  of  P  (e.g., 
for  r^  of  the  order  of  100  A,  P  =  10).  The  bold  assertion  made  in  applying  this  pore 
structure  model  to  soil  is  that  the  pore  size  distribution  is  valid  between  r^  and  r„,33j. 

To  investigate  the  role  of  the  measurement  scale  size  on  permeability,  consider  the 
largest  pore  r„,g^  contained  in  the  spherical  sample  of  radius  "a"  as  given  by  Eq.  (1).  Since 
in  Eq.  (31)  for  the  permeability  represents  the  largest  pore  in  the  medium,  the 
corresponding  value  of  "a"  is  denoted  a,„„  and  represents  the  largest  sample  size  for  which 
the  pore  sizes  will  scale  with  the  dimensions  of  the  sample.  From  Eqs.  (1)  and  (31) 


the  following  physical  characteristics: 

V  =2  cm/hr 
k  =0.6  Darcy 
0  =50% 

K<,=5 
p  =12 


04/3  V 


(32) 


Each  sphere  of  radius  a,„,^  will  contain  one  pore  of  size  A  20  x  20  grid  of  these  spheres 
will  be  characterized  by  the  dimension  40a„,a,t  and  contain  400  pores  of  size  r^.  Each  of 
these  pores  possess  probability  P,(rp)  of  being  interconnected.  Following  Eq.  (23),  P,(rp)  is 
approximately  0.0025  for  rp  sufficiently  close  to  r^.  Hence,  only  one  of  the  400  largest 
pores  in  this  20  x  20  grid  will  be  interconnected  and  the  error  in  the  measurement  of  the 
permeability  will  correspond  to  the  statistical  error  of  100%  associated  with  that  of  a  sample 
number  of  unity.  Carrying  this  argument  to  a  200  x  200  grid  of  dimension  400a,nax  >  there  will 
be  100  interconnected  pores  corresponding  to  a  statistical  error  of  10%.  Similarly,  a  grid  of 
scale  40(X)an,a,5  will  reduce  the  error  to  1  % . 


Figure  3  illustrates  the  predicted  permeability  measurement  error  associated  with  the 
soil  sample  characterized  above  =0.3  cm).  Note  that  the  errors  associated  with  the 
measurement  of  permeability  become  negligible  as  the  measurement  scale  size  approaches 
several  meters.  This  has  been  confirmed  by  the  infiltration  data  of  Shouse  et.  al.  (1994).  The 
measured  value  of  hydraulic  conductivity  asymptotes  to  2  cm/hr  at  measurement  scales 
greater  than  4  meters.  At  smaller  measurement  scales,  the  inferred  measurement  error  is 
calculated  under  the  assumption  that  the  asymptote  is  precisely  2  cm/hr.  The  excellent 
agreement  between  the  predicted  and  inferred  error  supports  the  extension  of  the  pore  tree 
model  to  describe  porous  permeable  media. 
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V.  BULK  GASEOUS  DIFFUSION  IN  PARTIALLY  SATURATED  MEDIA 


The  extended  pore  tree  model  is  readily  adapted  to  partially  saturated  media  through 
the  assumption  that  all  of  the  water  is  contained  in  pore  sizes  between  r^j^  and  r,a,  while  only 
gas  is  contained  between  r„,  and  r^.  Since,  by  Eq.  (6),  porosity  is  distributed  as  In(rp)  in 
pore  space,  the  porosity  associated  with  the  air  filled  pores  (6J  is  approximated  by 


ft  —  ^  max*  sat'  p 


(33) 


where  r,^  is  treated  as  an  independent  variable  of  the  saturated  pore  structure.  No  gaseous 
diffusion  is  allowed  within  rp  ^  r,.,.  It  is  demonstrated  that  the  gas  diffusivity  scales  as  l/r,„ 
and  it  is  the  sensitivity  of  r^g,  to  the  saturated  volume  that  controls  the  saturated  diffusivity. 

The  diffusive  mass  flux  in  a  single  pore  is  given  by 


g  p 


-^Pc 

dx 


(34) 


where  Dg  is  the  continuum  gas  diffusion  coefficient  (Dg=0.2  cm^/s).  Just  as  in  the  case  of 
convection,  it  must  be  determined  whether  the  mass  flux  across  plane  A  A  in  Fig.  2  is 
determined  by  the  interconnectivity  of  the  smaller  pores  out  of  the  plane  (Fig.  2b)  or  by  only 
those  pores  that  are  interconnected  in  the  plane  (Fig.  2a).  Consider  any  pore  of  radius  r,  in 
plane  AA  of  Fig.  2b  to  be  the  trunk  of  a  tree.  Each  pore  of  size  rp  within  the  tree  possesses 
the  probability  P/(r^)  of  being  interconnected  and  each  interconnected  pore  in  the  tree  will 

carry  the  mass  flow  rate  .  Since  there  are  N,drp  (Eq.  10:  N,=r,Vrp‘‘)  pores  in  size 

range  drp  within  the  tree,  the  total  mass  flow  rate  through  each  and  every  trunk  of 

radius  r,  in  plane  AA  becomes 


N,  dr^ 


(35) 


Integration  of  Eq.  (35)  yields  the  mass  flux  (case  2b)  for  each  trunk  of  radius  r. 


I  ; 


(36) 


where  a  ln(rggg,,/rgj  term  was  eliminated  via  Eq.  (33)  and  it  has  been  assumed  that  r,  >  >  r,,,. 
This  introduces  an  error  for  the  smaller  trees  close  to  r,,,  but  since  bulk  diffusion  (Db^j^ )  will 
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be  shown  to  be  dominated  by  the  largest  trees,  the  approximation  is  valid. 

If  the  mass  flux  through  plane  AA  is  limited  by  the  pores  that  are  interconnected  in 

that  plane  (case  2a),  the  mass  flux  is  expressed  as  Pjif’s) ,  and  it  is  immediately  seen 

that 


MAr^)>Mpir^PM 


(37) 


i.e.,  case  2b  dominates  case  2a.  Since  the  mass  diffusion  in  plane  AA  is  determined  by  the 
pore  interconnectivity  of  the  smaller  pores  outside  of  plane  AA,  the  saturation  of  those 
smaller  pores  becomes  an  important  element  in  the  bulk  gaseous  diffusion. 

Since  each  pore  of  radius  r,  in  plane  AA  carries  mass  flux  ,  the  bulk  diffusion 

coefficient  is  obtained  by  integrating  Eq.  (36)  over  all  pores  in  that  plane 


"""  %Der^  _ 


(38) 


subject  to  the  approximation  that  the  gas  is  contained  in  the  largest  pores  (r,^  >  >  r„t ). 


D 


buJk 


D„QQ^r„„ 

g  a  max 


(39) 


Eliminating  r^  via  Eq.  (31),  the  bulk  diffusivity  Dbuu,  is  expressed  in  terms  of  the 
permeability  k. 


D 


bulk 


(40) 


The  bulk  diffusivity  cannot  increase  indefinitely  with  increasing  permeability  as 
inferred  by  Eq.  (40).  In  deriving  this  expression,  the  mass  flux  through  the  interconnected 
branches  of  the  tree  was  not  constrained  firom  exceeding  the  diffusive  capabilities  of  the  trunk 
itself.  To  correct  this  potential  problem,  the  limit  of  is  determined  as  the  maximum 
diffusive  flux  (i.e.,  Dg)  in  pore  r^  integrated  over  all  pores  in  plane  AA. 


^.0 


(41) 
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D  is  indicative  of  existing  models  for  the  diffusivity  within  a  completely  connected 
pore  structure  and  is  illustrated  in  Figure  4  together  with  the  predicted  values  of  Dbuik  for  an 
extended  range  of  values  of  permeability  and  the  saturation  radius,  r,^. 

Model  predictions  correspond  to  the  measured  values  of  0a=O.2  and  6=0.5  from 
Washington  et.  al.,  (1994),  and  the  diffusivity  data  suggest  a  value  of  r^  in  the  range  of 
10  pm  to  100  pm.  An  exact  comparison  of  the  present  theory  to  the  least  squared  fit  obtained 
by  Washington  et.  al.,  (1994),  suggests  a  value  of  30  pm.  This  least  squared  fit  also 
demonstrates  a  permeability  dependence  of  the  power  0.53  whereas  the  model  predicts  0.5. 

While  the  excellent  agreement  of  the  model  with  the  data  of  Washington  et.  al.,  (1994) 
does  substantiate  the  present  theory,  it  is  apparent  that  there  is  a  very  wide  range  of  possible 
values  for  Dbun,  in  partially  saturated  soils  which  will  depend  upon  an  unknown  saturation 
radius.  A  two  order  of  magnitude  decrease  in  the  saturation  radius  will  increase  the  bulk 
diffusivity  by  two  orders  of  magnitude  and  yet  the  corresponding  increase  in  the  air  filled 
porosity  is,  by  Eq.  (33),  only  33%.  Hence,  field  measurements  of  the  unsaturated  volume  are 
not  sufficiently  accurate  to  correlate  bulk  diffusivities.  If  bulk  diffiisivities  are  to  be  correlated 
with  field  data,  such  measurements  should  attempt  to  measure  the  saturation  radius. 
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VI.  HYDRODYNAMIC  DISPERSION 


The  distribution  of  velocity  (or  volume  flux)  with  pore  size  has  been  utilized  to 
develop  a  description  of  hydrodynamic  dispersion  (Simons,  1996b).  The  volume  flux 

through  a  single  pore  of  radius  r.  subjected  to  a  pressure  gradient  —  is  given  by  Eq.  (24). 

dx 


8|i  1 


The  corresponding  fluid  velocity  v^(rp  within  that  single  pore  is 


(42) 


(43) 


and  it  is  evident  that  the  fluid  within  the  larger  pores  will  convect  ahead  of  the  fluid  in  the 
smaller  pores  inducing  an  apparent  difiusional  process  commonly  referred  to  as  hydrodynamic 
dispersion. 

It  has  been  shown  (Eq.  28)  that  the  volume  flux  QJr^  within  all  pores  of  radius  rp 
across  a  fixed  plane  is  limited  by  the  interconnectivity  P/zp  of  the  pore  structure 


<?.(»>)  >  <?,(/>)  P/rp 


(44) 


and  the  total  volume  flow  rate  Q  across  the  cross  sectional  area  A  of  the  pore  structure  is 
obtained  by  integrating  Eq.  (44)  over  all  pores  in  that  cross  section.  Using  the  pore  size 

distribution  function  denoted  by  g(r^ ,  Q  becomes 


^jnioi 

from  which  permeability  k  is  defined  as 


(45) 


A{-dpldx) 


6^ 
16p  , 


(46) 
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and  the  mean  pore  velocity  v  is  defined  as  the  average  fluid  velocity  over  the  porous  area. 


^-7-  Q  _ki-dpldx)  (47) 

QA  0p 

The  mean  velocity  associated  with  the  pore  radius  Tp  is  denoted  by  v^(rp  and  it  is  the 
variation  of  v^(rp  about  v  that  causes  the  "diffusion"  or  hydrodynamic  dispersion  of  a 

specific  species  about  the  mean  convection  velocity  v  .  The  aggregate  fluid  velocity  from 
plane  to  plane  must  be  determined  as  a  function  of  pore  radius.  In  any  given  cross  sectional 
plane,  the  fraction  Pj[r^  of  the  pores  were  determined  to  be  interconnected  and  carry  fluid 
velocity  v^(rp  while  the  non-interconnected  pores  in  that  size  range  do  not  support  a 
significant  fluid  velocity.  This  was  interpreted  as  an  effective  fluid  velocity  v^rp  P/rp 
occurring  in  all  pores  of  radius  rp  in  a  given  cross  sectional  plane.  An  expression  for  v^(rp 
is  utilized  which  scales  with  v^rp  P/rp  and  recovers  the  mean  convection  velocity  v  in  the 
limit  of  f‘p=^p,  where  is  the  pore  size  associated  with  velocity  v  and  is  yet  to  be 
determine. 


(48) 

^  Mr^lrp) 

The  mean  velocity  associated  with  pores  of  size  rp  will  be  used  to  relate  a  species 
concentration  profile  in  pore  space  to  that  in  physical  space  using 

^  (49) 

V  V  f 

where  x  is  a  spatial  coordinate  in  the  concentration  profile.  That  profile  will  contain  the 
concentration  c  moving  with  mean  velocity  v  and  associated  with  pore  size  were  both 
c  and  are  yet  to  be  determined.  The  system  will  lead  to  the  "diffusion"  or  hydrodynamic 
dispersion  of  a  specific  species  about  the  concentration  c  moving  with  the  mean  convection 
velocity  v  . 

From  Eq.  (48)  it  is  seen  that  the  mean  velocity  associated  with  pores  of  size  rp 
possesses  a  relative  maximum  at  r„„l^  as  illustrated  in  Figure  5.  Let  r,  and  T2  represent 
the  two  values  of  rp  corresponding  to  a  particular  value  of  v^(rp .  Clearly,  the  fluid  in  the 
pore  of  radius  convects  ahead  of  that  in  the  pores  of  radius  rj  and  rj.  Consider  flow 
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in  a  saturated  porous  medium  in  which  all  fluid  behind  a  fixed  plane  normal  to  the  flow  is 
suddenly  injected  with  red  dye.  All  dyed  fluid  in  pores  whose  radius  is  between  ri*  and  r2* 
will  translate  with  velocity  equal  to  or  greater  than  v„*  and  move  to  a  corresponding  location 
X*  or  greater.  The  concentration  profile  c(x)/Co  of  the  red  dye  at  location  x  will  scale  with  the 
volume  flux  whose  velocity  is  greater  than  v„(x)  from  within  pores  whose  radius  is  between 
ri(Vm(x))  and  rjfv^fx)).  From  these  arguments,  it  follows  that  the  normalized  concentration 
profile  is  given  by 


=  / <?.(V 

'■l  'll* 


(50) 


where  ri(v„(x))  and  rjfv^fx))  are  related  through  Eq.  (48), 


(51) 


and  c(x)/Co  becomes 


cix)lc^  =  {rl-rl)lrl^ 


(52) 


The  concentration  profile  of  the  red  dye  in  }X)re  space  is  readily  obtained  by  solving  Eq.  (51) 
for  and  rj/r^^^  parametric  in  terms  of  r/rj  between  zero  and  unity.  The  corresponding 

dye  concentration  profile  in  pore  space  is  illustrated  in  Figure  12.  Note  that  the  pore 
containing  the  largest  velocity  corresponds  to  the  "front"  of  the  spatial  profile  (Figure  6) 
where  c=0  and  the  pores  containing  the  lower  velocities  are  at  the  rear  of  the  spatial  profile 
where  c  tends  to  unity. 
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Figure  12.  Concentration  in  Pore  Space 


Figure  13.  Hydrodynamic  Dispersivity^ 
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Eq.  (50)  represents  the  concentration  profile  of  the  red  dye  in  pore  size  space,  pore 
velocity  space,  and  physical  space.  The  concentration  c  associated  with  the  mean  velocity  v 
may  be  determined  from  the  conservation  of  the  flux  of  the  concentration  profile.  This  is 
expressed  as 


/  c] 


0.(/>)  gir^  A  dr^ 


(53) 


where  the  integral  is  split  into  integrals  over  the  ri  and  T2  branches  and  obtained  numerically 
firom  the  profile  illustrated  in  Figure  12.  Integration  yields 


clc^=0.5  (54) 

It  is  seen  from  Figure  12  that  there  are  two  pore  radii  in  pore  space  that  satisfy 
Eq.  (54).  On  the  ri  branch  ^  is  0.38332  r„„„  and  on  the  rj  branch  we  obtain  0.80428  r^nj,- 
In  each  case 


ln(r^/rp  =  0.1409 


(55) 


and  Eqs,  (48)  and  (49)  assign  a  unique  value  of  x  to  each  value  of  rj  and  rj.  We  define 


^  rl  ln(r^/ri) 

vt  2  0,1409 

'max 


(56) 


and/or  the  identical  expression  in  r2.  The  concentration  profile  may  be  expressed  in  physical 
space  as  a  function  of  q  corresponding  to  that  in  pore  space  as  a  function  of  rj  and  rj.  The 
concentration  profile  in  physical  space  is  illustrated  in  Figure  6  where  all  red  dye  would  be 
located  at  the  location  q  =  1  were  it  not  for  the  velocity  variation  within  pore  size.  The 
"front"  of  the  concentration  profile  has  progressed  30%  ahead  of  the  mean  while  the  slower 
fluid  in  the  smaller  pores  has  trailed  considerable. 

If  we  were  standing  at  the  end  of  a  porous  flow  tube,  the  concentration  passing  our 
location  would  exhibit  a  time  dependence  (  c  vs.  1/q)  illustrated  in  Figure  7.  The  faster  fluid 
breaks  through  ahead  of  the  mean  flow  (c/c^=0.5)  while  the  fluid  in  the  smaller  pores  is 
delayed.  Data  of  Elrich,  et  al.,  (1966),  as  reported  by  Brusseau  and  Rao  (1989),  are 
illustrated  in  Fig.  8  and  support  the  current  theory  of  hydrodynamic  dispersion.  While  the 
long  "tail"  of  the  concentration  profile  is  generally  attributed  to  nonequilibrium,  the  current 
pore  structure  model  attributes  it  to  the  wide  pore  size  distribution  occurring  in  soil.  The 
narrow  pore  size  distributions  occurring  in  laboratory  "soil"  will  yield  a  much  narrower 
concentration  profile,  i.e.,  less  hydrodynamic  dispersion. 
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The  local  value  of  the  dispersivity  required  for  the  concentration  profile  to  be 
compatible  with  the  diffusion  equation  reflects  the  virtue  of  expressing  hydrodynamic 
dispersion  as  a  Fickian  process.  The  hydrodynamic  dispersion  (D^  in  length^/time)  normalized 

by  the  mean  velocity  v  is  the  dispersivity  X.  The  dispersivity  normalized  by  the  mean 
convection  distance  v  f  is  expressed  in  terms  of  the  concentration  profile  via 

X  _  ( 1  ~'n )  dcjdi] 
vt  ^cjdr\^ 

and  illustrated  in  Figure  13.  Figure  13  illustrates  that  neither  X  nor  Xjvt  is  constant,  i.e., 
the  concentration  profile  could  never  be  described  by  a  Fickian  process.  Future  research  will 
be  directed  toward  using  this  analytic  pore  structure/  pore  transport  model  to  help  derive  an 
analytic  expression  for  hydrodynamic  dispersion  that  replaces  the  Fickian  process  in  the 
transport  equations. 
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vn.  COUPLING  SUBSCALE  REACTIONS  TO  MACROSCOPIC  TRANSPORT 


Analysis  of  the  permeability,  bulk  diffusivity,  and  hydrodynamic  dispersion  has 
utilized  the  interconnectivity  of  the  pores  to  determine  the  distribution  of  the  convection 
velocity  with  pore  size.  This  analysis  suggests  the  existence  of  a  permeable  sub-range  in  the 
pore  structure  which  does  not  contribute  significantly  to  the  bulk  permeability  but  in  which 
convection  dominates  diffusion.  It  is  the  balance  of  the  sub-scale  convection  with  the  small 
scale  diffusion  that  will  control  contaminant  transport  and  in-situ  remediation.  The  size 
distribution  of  the  pores  and  grains,  and  the  variations  in  fluid  velocity  within  and  between 
pores  of  different  sizes  is  critical  to  interfacing  the  transport  processes.  A  methodology  is 
developed  to  couple  subscale  diffusion,  convection  and  chemical  reactions  to  the  macroscopic 
transport  in  order  to  accurately  describe  contaminant  transport  and  in-situ  remediation  in 
Ground  Water  Simulation  codes. 

A  typical  species  transport  equation  in  permeable  soil  is  of  the  form: 


dc  _  2)  +  cM 

dt  dx  ^  dx^  M 


(58) 


where  c  is  the  local  species  concentration,  v  is  the  mean  convection  velocity  as  determined 
by  the  local  pressure  gradient  and  permeability  k,  and  is  the  hydrodynamic  dispersion. 
The  dc/dt  term  is  the  true  unsteady  term  and  the  source  and/or  sink  of  species  c  due  to 
chemical  and/or  physical  processes  within  the  subscale  pore  structure  is  written  as  cAf/Af 

where  MJM  represents  a  bulk  rate  (1/t)  of  production  or  consumption.  The  bulk  rate  Af/Af 
will  be  related  to  gradients  of  c  on  length  scales  of  order  millimeters  within  the  subscale  pore 
structure  and  cannot  be  expressed  in  terms  of  the  dc/dx  of  the  macroscopic  transport  grid. 

To  describe  the  subscale  chemical  reactions  and  transport,  a  "grain”  of  soil  of  radius 
2Lg  is  isolated  from  the  rest  of  the  medium.  The  grain  size  will  be  chosen  sufficiently  small 
that  molecular  diffusion  and  coupled  chemical  reactions  will  dominate  the  subscale  transport. 
The  pore  tree  was  developed  (Simons,  1982,  1983a,  1988)  to  treat  coupled  chemical  reactions 
and  molecular  diffusion  within  immobile  porous  grains.  A  reactant  species  will  diffuse  into 
the  pores  and  react  on  the  walls  of  the  pore.  The  chemical  reaction  on  the  walls  of  the  pore 
establishes  a  species  gradient  which  is  greater  than  c/ag  and  accelerates  the  molecular 
diffusion.  The  effective  molecular  diffusion  coefficient  is  the  square  root  of  the  product  of  the 

molecular  diffusivity  and  the  chemical  rate  constant.  The  solution  for  M  in  this  limit  scales  as 
a^^  because  transport  is  limited  by  diffusion  across  the  surface  area  ag^  while  the  species 
gradient  is  independent  of  ag.  The  particular  form  of  the  solution  for  contaminant  transport 
and  soil  remediation  will  vary  with  the  number  of  reactions  and  kinetic  mechanisms  but  the 

fundamentals  of  the  problem  depict  that  M  will  scale  as  ag^  in  this  diffusive  limit. 

As  the  grain  size  is  reduced,  the  diffusive  abilities  of  the  grain  are  able  to  keep  up 
with  the  kinetics  of  the  chemical  reactions.  The  species  concentration  gradient  becomes  small 


25 


compared  to  c/ag  and  the  species  reaction  rate  scales  with  the  volume  of  the  reactant.  In  this 

"kinetic  limit,"  ik  will  scale  as  a.^.  Since  the  local  mass  of  the  reactant  species  M  also 

scales  as  ag^  it  follows  that  M/M  scales  as  ag®  in  the  kinetic  limit  and  as  ag  ‘  in  the 
diffusive  limit.  This  behavior  is  expressed  as 


—  =  .Rj  ( Kinetic  Limit  :ag<R^/  R^) 


and 


(59) 


JVm.  *  *  * 

—  =J^/a^  (Diffusive Limit :ag>R2 1  Ri) 


(60) 


where  it  is  assumed  that  R^  and  R^  may  be  rigorously  determined  for  specific  kinetic 
processes  in  soil  through  direct  application  of  the  pore  tree/pore  transport  model. 

The  reactant  species  must  be  supplied  to  the  edge  of  the  grain  at  a  rate  compatible 
with  the  rate  of  species  consumption  within  the  grain.  Convection  in  the  subscale  will 
transport  the  species  across  the  cross  section  Tiag^  at  a  rate  depicted  by  the  permeability  kg 

which  is  appropriate  to  the  length  scale  ag.  Since  M  is  and  Q  may  be  obtained 

from  equation  (25),  M/M  becomes 

M  ^  c  p(kg/\i)(nag)(-dpldx) 

^  cpe(4%ag/3) 

where  p  is  the  density  of  the  reactant  species  (gas  or  liquid  phase).  The  integral  for  the 
permeability  (Eq.  30)  is  dominated  by  the  largest  pores.  If  r„,g  is  the  largest  pore  in  the  grain 
and  r„„„j  is  the  largest  pore  in  the  soil,  kg  appropriate  to  the  grain  will  scale  with  k,  the 
permeability  of  the  soil,  via 


k 

'max 

where  it  is  assumed  that  r„,„/  r,„g  is  roughly  half  the  width  of  the  entire  pore  distribution.  i.e. , 
the  (1  +21n)  term  is  approximately  equal  to  3  (Eq.  4). 

The  radius  (r„,g)  of  largest  pore  in  the  grain  of  radius  Ug  is  directly  proportional  to  ag 
just  as  r^  is  directly  proportional  to  an,,*,  the  radius  of  the  largest  grain  in  which  the  pore 
size  will  scale  with  the  dimension  of  the  sample.  Hence,  equation  (62)  becomes 
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(63) 


^  P 

and,  using  Eq.  (32)  for  a„,ax  >  E<1-  (63)  is  rewritten  as 

*  S16  Kl^k 

such  that  the  convective  limit  of  M/M  follows. 


(64) 


M  ^  0 — (  dpidx)  ^  a  {Convective limit)  (65) 


The  convective  limit  of  M/M  does  not  explicitly  depend  on  the  bulk  permeability  k  because 
k  is  controlled  by  the  largest  pores  in  the  soil  and  cannot  distribute  the  reactants  on  a  scale  as 
fine  as  these  small  grains. 

These  results  are  illustrated  in  Figure  9  and  indicate  that  the  convective  limit  increases 
with  larger  grains  while  the  diffusive  limit  decreases  with  increasing  grain  size.  The 
convective  and  diffusive  limits  are  equal  at  ag  =  a^q 

where  it  is  apparent  that  for  grains  smaller  than  a^q  ,  the  convection  cannot  keep  up  with  the 
diffusion  and  for  grains  larger  than  a^q,  the  diffusion  cannot  keep  up  with  the  convection.  The 

system  will  attain  a  natural  equilibrium  at  ag  =  and  the  corresponding  value  of  M/M 
becomes 


-  Smaller  Of  (67) 

JM. 

where  it  has  been  noted  that  the  solution  at  ag  =  a^q  cannot  exceed  the  limit  of  kinetic  control. 

The  above  expression  for  the  bulk  source  and/or  sink  of  species  c  due  to  chemical 
and/or  physical  processes  within  the  subscale  pore  structure  is  incomplete  without  specific 

models  for  and  .  Future  research  will  be  directed  toward  developing  a  library  of  such 
models.  The  most  basic  processes  of  contaminant  diffusion  and  adsorption  are  briefly 
described  below. 
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Vin.  IMMOBILE  SOLUTIONS:  TWO  LIMITING  CASES 


Vni.l  General  Approach 

The  pore  structure  model  described  above  has  been  coupled  to  a  transport  model 
(Simons,  1982)  that  is  capable  of  describing  the  simultaneous  action  of  diffusive  and 
kinetically  controlled  processes  in  the  immobile  region.  The  pore  transport  model  was 
developed  to  describe  quasi-steady  char  oxidation  in  which  both  adsorption  and  desorption 
occur  simultaneously.  This  transport  model  will  be  adapted  to  describe  the  transient 
chemisorption  and  remediation  processes  in  soil  such  that  we  will  be  able  to  couple  the 
subscale  reactions  to  the  macroscopic  transport  as  described  above. 

Each  pore  that  reaches  the  exterior  surface  of  the  particle  is  depicted  as  the  trunk  of  a 
tree.  The  size  distribution  of  tree  trunks  on  the  exterior  surface  of  the  particle  is  denoted  by 

g{r^%a}drj  where  is  functionally  identical  to  .  Each  trunk  of  radius  rj  is 
associated  with  a  specific  tree-like  structure  as  illustrated  in  Fig.  1,  A  pore  structure  with 
continuous  branching  to  ever  decreasing  pore  radii  is  depicted  as  being  attached  to  a  uniform 
trunk  of  radius  r,.  The  radius  and  number  of  pores  is  a  unique  function  of  the  distance  x  into 
the  tree.  The  coordinate  x  is  skewed  in  that  it  follows  a  tortuous  path  through  the  branches  of 
the  tree.  Let  n(x)  represent  the  number  of  pores  of  radius  rp  at  location  x  in  a  tree  of  trunk 
radius  r^.  An  analysis  (Simons,  1982)  of  this  pore  tree  has  demonstrated  that 


n(.x)=rilTl{x) 

and  the  coordinate  x  is  related  to  tp  and  the  length  of  the  tree  trunk  by 


(68) 


drjdx=-rjl^ 


(69) 


When  this  porous  structure  is  placed  in  an  environment  with  a  reactive  gas,  the 
reactive  gas  will  diffuse  into  the  pore  tree  and  react  with  the  material  that  constitutes  the 
walls  of  the  pore.  The  diffusion  of  the  reactive  gas  through  n  pores  of  radius  Tp  is  balanced  by 
the  continuous  reaction  of  the  gas  at  the  walls  of  the  pore.  We  write 


where  c  is  the  mass  fraction  of  the  reactant  gas,  Dg  is  the  self  diffusion  coefficient  of  the 
reactant  species,  Pg  is  the  combined  density  of  all  gases  in  the  pore,  Pg  is  the  combined 
pressure  (in  atmospheres)  of  all  gases  in  the  pore,  and  k,  is  the  rate  constant  (g/s/cmVatm)  for 
a  continuous  first  order  reaction  which  is  a  ftinction  of  particle  temperature  alone. 
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The  total  reaction  rate  of  the  pore  tree,  ,  is  related  to  the  gradient  of  c  at  x=0  by 


dx 


(71) 


However,  to  obtain  the  value  of  dc/dx  at  x=0,  Eq.  (70)  must  be  integrated  subject  to  the 
boundary  conditions  that  c=Co  at  x=0  (rp  =r, )  and  dc/dx=0  at  x=Xt  (rp  =r,„J.  The  total 

reaction  rate  of  the  porous  particle  is  obtained  by  integrating  over  all  trees 


Mj,=  f  M^4%a^g(r)dr^  (72) 


The  pore  tree  model  has  been  used  to  successfully  describe  char  oxidation 
(Simons,  1979b;  Lewis  and  Simons,  1979;  Simons,  1982  &  1983a),  coal  pyrolysis  (Simons, 
1983b  &  1984)  and  the  catalytic  cracking  of  benzene  by  porous  iron  oxides  (Simons  et  al., 
1986).  It  was  also  used  to  successfully  describe  sulfur  sorption  (SO2  and  H2S)  by  porous 
calcine  (CaO)  in  the  limit  of  zero  utilization  (Simons  and  Rawlins,  1980;  Simons  et  al.,  1984) 
and  was  later  extended  to  include  CaS04  and  CaS  deposits  (Simons  and  Garman,  1986; 
Simons  et  al.,  1987;  Simons,  1988).  The  subsequent  determination  of  the  controlling  physical 
parameters  led  to  a  new  concept  for  the  optimization  of  the  sulfur  sorption  process  (Simons, 
1991;  Simons  et  al.,  1992)  through  spray  drying  of  water  soluble  organic  calcium  solutions  to 
control  the  sorbent  pore  structure.  This  transport  model  will  be  adapted  to  describe  the 
transient  chemisorption  and  remediation  processes  in  soil  such  that  we  will  be  able  to  couple 
the  subscale  reactions  to  the  macroscopic  transport  as  described  above. 


VII1.2  Contaminant  Adsorption  in  Saturated  Media 

Consider  a  contaminant  in  a  completely  saturated  medium.  The  contaminant  will 
diffuse  into  the  pore  structure  (Figure  1)  as  described  by  Eq.  (70).  Assume  that  the  diffusion 
is  rate  limiting  and  the  kinetic  term  (adsorption  rate)  in  Eq.  (70)  is  relatively  fest.  Within  this 
idealization,  the  walls  of  the  pores  between  x=0  and  x=X(.(t)  will  be  completely  contaminated 
while  the  walls  of  the  pores  for  x  >  x<.(t)  will  not  be  contaminated.  All  of  the  reaction  is 
assumed  to  occur  at  (at  and  around)  x<,(t)  where  Xj,(t)  expands  throughout  the  entire  pore  tree 
as  time  approaches  infinity.  The  first  integral  of  Eq.  (70)  between  x=0  and  x=X(.(t)  becomes 

=  (73) 

where  denotes  the  liquid  density,  denotes  the  contaminant  diffusion  coefficient  and 
denotes  the  net  mass  flux  of  the  contaminant  into  the  pore  tree  as  a  function  of  time. 
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Using  n(x)  and  dfp/dx  defined  above,  Eq.  (73)  may  be  integrated  subject  to  the 
boundary  condition  that  c=Co  at  x=0  (rp=rd  and  c=0  at  x=Xc  (rp=r<,).  The  contaminant 
concentration  in  the  liquid  within  the  pore  structure  becomes 


Hr.lO 


and  the  net  mass  flux  of  the  contaminant  into  the  pore  tree  is  given  by 


(74) 


MM  = 


(75) 


where  the  ln(rp/r^)  term  has  been  replaced  with  ^  (Eq.  4). 

The  behavior  of  M^(x^)  near  complete  contaminant  adsorption  (for  r^  near  r,^  is 
incorrect  and  an  exact  solution  to  Eq.  (70)  is  required  to  precisely  simulate  this  limit. 
However,  Eq.  (75)  may  be  modified  to  more  closely  simulate  the  mass  flux  in  the  limit  of 
complete  contaminant  adsorption.  If  the  walls  of  the  pore  tree  adsorb  contaminant  at  the  rate 
depicted  by  Eq.  (75),  the  reaction  will  proceed  until  the  entire  surface  area  S,  of  the  pore  tree 
is  coated  with  contaminant.  Denoting  the  maximum  surface  contaminant  adsorption  level  as 

(mass  per  unit  area),  the  characteristic  time  x  over  which  Eq.  (75)  is  valid  is  SJM,. 
Using  Eq.  (13)  for  the  surface  area  Sj  of  the  pore  tree,  x  becomes 

t  .  (76) 

C  Dr  Dr  r  . 
o  L mm 


and  is  rewritten  as 


MM 


/.P 


exp(-tjx) 


(77) 


The  total  adsorption  rate  of  the  porous  grain  is  obtained  by  integrating  M,  over  all 
trees  as  in  Eq.  (72)  above.  At  late  time,  the  integral  is  dominated  by  the  largest  pore  trees  in 
the  grain.  The  radius  of  the  largest  tree  trunk  is  related  to  the  grain  size  by  equation  (1) 
and  the  total  adsorption  rate  of  the  porous  grain  is  see  to  be 


3  It  0  Og 

P^ 


(78) 
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As  contaminant  adsorption  nears  completion,  the  total  mass  adsorbed  approaches  the 
product  of  the  maximum  surface  contaminant  adsorption  level  (mass  per  unit  area)  and  the 
total  internal  surface  area.  Hence 


3pr_ 


(79) 


and  the  immobile  solution  for  M/M  becomes 


M  _  ^  ^I.  P  T.  _  ^4 

"  8  p  oy.  ^ 


(80) 


where 


^  9c^DiPi,r^ 


(81) 


and  the  coupled  mobile/immobile  solution  is  developed  as  described  in  Section  Vn 


(82) 


where  the  viscosity  p  in  1^  reflects  the  properties  of  the  liquid  in  the  saturated  medium. 


One  feature  of  the  immobile  solution  is  immediately  obvious:  M/  M  scales  as  l/a^^ 
instead  of  l/a^  as  determined  for  the  coupled  reaction/diffusive  systems  appropriate  to 
combustion  applications.  This  occurs  because  the  late  time  concentration  gradient  within  the 
grain  is  and  introduces  an  additional  1/ag  to  the  scaling  relation.  It  is  emphasized  that 
Eq.  (80)  is  valid  only  at  late  time  as  the  contaminant  adsorption  reaches  the  total  mass 

adsorption  limit.  At  early  time  the  integral  for  M  is  dominated  by  the  smallest  trees,  the 

corresponding  concentration  gradient  is  much  greater  than  and  M  is  r„,g/r„,i„  greater  than 
Eq.  (80).  The  early  time  solution  then  scales  as  l/Ug  but  represents  the  adsorption  rate  only  at 
extremely  low  fractions  of  the  total  mass  adsorption  limit.  This  clearly  indicates  that  a 
complete  immobile  model  for  this  process  will  involve  a  numerical  subroutine  rather  than  a 
simple  analytic  expression.  While  such  a  model  is  quite  tractable  and  would  be  coupled  to  the 
mobile  solution  as  previously  described,  it  is  beyond  the  scope  of  the  current  study. 


Vin.3  Contaminant  Adsorption  in  Unsaturated  Media 


Consider  a  contaminated  gas  in  a  completely  unsaturated  medium.  The  contaminant 
will  diffuse  into  and  deposit  onto  the  walls  of  the  pores  in  exactly  the  same  way  as  described 
for  the  saturated  case.  Within  the  approximation  that  the  kinetics  of  the  adsorption  are  fast 
with  respect  to  the  diffusion,  the  immobile  solution  is  given  by  Eq.  (80)  with  the  gas  density 
pg  replacing  the  liquid  density  Pl  and  the  gas  phase  contaminant  diffusion  coefficient  Dg 
replacing  the  liquid  phase  contaminant  diffusion  coefficient  .  The  immobile  solution  is 
written  as 


M  _  ^ 


M 


(83) 


where 


8Po, 


(84) 


and  the  coupled  mobile/immobile  solution  is 


(85) 


where  the  viscosity  u  in  now  reflects  the  properties  of  the  gas  in  the  unsaturated  medium. 

It  is  again  emphasized  that  Eq.  (85)  is  valid  only  at  late  time  as  the  contaminant 
adsorption  reaches  the  total  mass  adsorption  limit  and,  just  as  in  the  saturated  case,  a 
complete  model  for  this  time  dependent  process  will  involve  a  numerical  subroutine  rather 
than  a  simple  analytic  expression.  However,  the  analytic  form  of  the  solution  appropriate  to 
late  time  offers  a  relatively  simple  expression  to  couple  into  the  Ground  Water  Simulation 
codes  for  the  purpose  of  testing  this  method  of  coupling  the  mobile  and  immobile  regions. 


K.  SUMMARY 


The  pore  tree  model  has  been  extended  to  describe  the  permeable  pore  structure  which 
characterizes  the  subsurface  transport  of  gas  and  water  in  soil,  the  dispersion  of  contaminants, 
and  the  in-situ  remediation  of  contaminated  sites.  The  random  nature  of  the  pore  structure, 
which  formed  the  basis  of  the  statistical  derivation  of  the  pore  tree,  is  applied  to  porous  soil 
and  sand.  The  interconnectivity  of  the  pore  structure  is  obtained  via  a  statistical  determination 
of  the  "branches"  that  are  common  to  several  trees  to  allow  convection  and  bulk  diffusion 
through  the  large  scale  (mobile)  structure  in  addition  to  diffusion  and  coupled  chemical 
reactions  within  the  smaller  scale  (immobile)  structure.  The  statistical  analysis  reported  above 
has  determined  that  the  probability  of  pore  interconnectivity  extends  across  the  entire  pore 
size  range,  with  a  slight  increase  in  the  probability  accompanying  a  decreasing  pore  size. 
While  permeability  is  dominated  by  the  largest  pores,  it  is  also  important  to  establish  the  level 
of  convection  and  diffusion  that  is  occurring  at  the  intermediate  scales  in  order  to  accurately 
relate  large  scale  bulk  transport,  intermediate  scale  convection,  small  scale  diffusion  and 
coupled  chemical  reactions. 

The  permeability  across  a  given  plane  is  limited  by  the  largest  pores  that  are 
interconnected  in  that  plane.  The  statistical  analysis  has  determined  that  approximately  one 
quarter  of  one  percent  of  all  large  pores  are  interconnected.  This  establishes  a  very  coarse 
grid  for  the  permeability  which  leads  to  measurement  scale  size  errors.  The  extended  pore 
tree  model  has  successfully  explained  the  measurement  errors  in  the  permeability  of  soil  due 
to  the  measurement  scale  size  (Shouse,  et.al.,  1994)  which  has  indirectly  confirmed  the  low 
probability  of  the  interconnectivity. 

The  bulk  gaseous  diffiisivity  across  a  given  plane  is  shown  to  be  limited  by  the 
interconnectivity  of  the  smaller  branches  outside  of  that  plane.  These  small  pores  may  be 
saturated,  resulting  in  a  strong  dependence  of  the  diffusivity  on  the  radius  of  the  saturated 
pore.  A  comparison  of  the  present  theory  to  the  diffusivity  data  of  Washington  et.  al.,  (1994) 
suggests  a  saturation  radius  of  30  pm.  While  the  excellent  agreement  with  the  data  does 
substantiate  the  present  theory,  the  diffusivity  in  partially  saturated  soil  is  very  sensitive  to  an 
unknown  saturation  radius.  If  bulk  diffusivities  are  to  be  correlated  with  field  data,  such 
measurements  should  attempt  to  measure  the  saturation  radius. 

The  permeability  and  the  bulk  diffusivity  have  tested  two  extreme  limits  of  the  pore 
structure  and  pore  interconnectivity  concepts.  Permeability  is  limited  by  the  in  plane 
interconnectivity  (Fig.  2a)  and  bulk  gaseous  diffusion  is  limited  by  the  out  of  plane  (Fig.  2b) 
interconnectivity.  Permeability  is  limited  by  the  large  pore  interconnectivity  and  bulk  diffusion 
is  limited  by  the  interconnectivity  of  the  smaller  pores.  The  apparent  success  of  these 
concepts  over  a  very  broad  pore  size  range  suggests  that  the  extended  pore  tree  model  may  be 
used  to  develop  the  subscale  convection  necessary  to  couple  the  mobile  and  immobile  regions. 

Since  hydrodynamic  dispersion  is  generated  by  the  fluid  velocity  differences  between 
the  smaller,  convection  dominated  pores  (Garabedian  et  al.,  1991),  determination  of  the 
hydrodynamic  dispersion  in  a  homogeneous  field  has  been  used  (Simons,  1996b)  to  test  the 
pore  interconnectivity/pore  structure  concepts.  Data  of  Elrich,  et  al.,  (1966),  as  reported  by 
Brusseau  and  Rao  (1989),  support  the  current  theory  of  hydrodynamic  dispersion.  While  the 
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long  "tail"  of  the  concentration  profile  is  generally  attributed  to  nonequilibrium,  the  current 
pore  structure  model  attributes  it  to  the  wide  pore  size  distribution  occurring  in  soil.  The 
narrow  pore  size  distributions  occurring  in  laboratory  "soil"  will  yield  a  much  narrower 
concentration  profile,  i.e.,  less  hydrodynamic  dispersion.  Future  research  will  be  directed 
toward  1),  evaluating  the  hydrodynamic  dispersion  as  a  function  of  pore  size  range  in  the  soil 
sample  and  2),  using  this  analytic  pore  structure/  pore  transport  model  to  help  derive  an 
analytic  expression  for  hydrodynamic  dispersion  that  replaces  the  Fickian  process  in  the 
transport  equations. 

Analysis  of  the  permeability,  bulk  diffusivity,  and  hydrodynamic  dispersion  has 
utilized  the  interconnectivity  of  the  pores  to  determine  the  distribution  of  the  convection 
velocity  with  pore  size.  This  analysis  suggests  the  existence  of  a  permeable  sub-range  in  the 
pore  structure  which  does  not  contribute  significantly  to  the  bulk  permeability  but  in  which 
convection  dominates  diffusion.  It  is  the  balance  of  the  sub-scale  convection  with  the  small 
scale  diffusion  that  will  control  contaminant  transport  and  in-situ  remediation.  The  size 
distribution  of  the  pores  and  grains,  and  the  variations  in  fluid  velocity  within  and  between 
pores  of  different  sizes  is  critical  to  interfacing  the  transport  processes.  A  methodology  has 
been  developed  to  couple  subscale  diffusion,  convection  and  chemical  reactions  to  the 
macroscopic  transport  in  order  to  accurately  describe  contaminant  transport  and  in-situ 
remediation  in  Ground  Water  Simulation  codes.  The  fundamental  approach  to  this  application 
was  described  in  Section  Vn  and  two  specific  examples  aj^ropriate  to  limiting  cases  of 
subscale  transport  were  developed  in  Section  VDI.  Future  research  will  be  directed  toward 
developing  a  library  of  models  for  chemical/diffiisive  processes  in  the  immobile  region  such 
that  they  may  be  coupled  to  the  mobile  region  via  the  technique  derived  above.  Such  models 
will  provide  Ground  Water  Simulation  Codes  with  physically  realistic  submodels  for 
contaminant  transport  and  remediation  studies  that  can  be  used  to  scale  data  from  the 
laboratory  to  the  field  and  from  one  field  site  to  another. 
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XI.  NOMENCLATURE 


A 

Ae 

a 


^max 


c 


Dh 

D, 

Dl 

Dbuk 

f(rp) 

g(r^ 

G, 

/(rp 

k 


K 

Ko 

Ip 

It 


Mprp 

M«(rp 

M/Af 

N. 

n(x) 

P 


<? 


<?prp 


fc 

^mg 

^flat 

ri 


r2 


arbitrary  cross  sectional  area 

area  within  A  that  is  connected 

radius  of  sphere  in  a  porous  medium 

maximum  value  of  a  in  which  pore  structure  is  self  similar 

radius  of  arbitrary  "grain"  of  soil 

species  mass  fraction  (pc/pg) 

mean  species  concentration 

hydrodynamic  dispersion  coefficient 

continuum  gas  phase  diffusion  coefficient 

contaminant  di^sion  coefficient  in  the  liquid 

gas  diffusion  coefficient  in  partially  saturated  porous  media 

gas  diffusion  coefficient  in  fully  connected  porous  media 

number  of  pores  per  unit  volume  of  radius  rp 

number  of  pores  per  unit  area  of  radius  rp 

number  of  pores  in  tree  r,  and  in  plane  A  of  radius  rp 

common  branch  distribution  function  I  =  gPj 

permeability 

permeability  of  grain  ag 

constant  relating  the  pore  length  to  the  radius 

length  of  a  pore  of  radius  rp 

length  of  a  tree  trunk  of  radius  r, 

diffusive  mass  flux  in  pore  of  radius  rp 

diffusive  mass  flux  in  tree  rp  limited  by  interconnectivity 

mobile/immobile  rate  (l/t) 
number  of  pores  in  tree  of  radius  rp 
number  of  branches  at  location  x 
probability  of  interconnectivity 
gas  or  liquid  pressure 

convective  volume  flow  rate  across  cross  sectional  area  A 

convective  volume  flow  rate  in  pore  of  radius  rp 

convective  volume  flow  rate  in  tree  rp  limited  by  interconnectivity 

mobile/immobile  rate  constants 

radius  of  a  pore 

radius  of  trunk  of  tree  r, 

radius  of  penetration  of  contaminant  adsorption 

maximum  pore  radius  in  grain  ag 

radius  of  largest  pore  that  is  saturated 

largest  pore  corresponding  to  v^(rj) 

smallest  pore  corresponding  to  Vg,(rj) 
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Sp  specific  internal  surface  area  (m^/g) 

S,(rt)  surface  area  of  pore  tree  with  trunk  radius  rt 

t  time 

V,(r,)  volume  of  pore  tree  with  trunk  radius  r, 

v^(rp  velocity  associated  with 

v^(rp  mean  value  of  v^(rp  in  interconnected  pore  space 

V  mean  velocity  across  area  0A 

X  axial  distance  in  pore,  spatial  coordinate 

x^  axial  penetration  of  contaminant  adsorption 

Greek 

p  In  (  r„,^  /  r„i„ ) 

Tj  dimensionless  coordinate:  Eq  (56) 

9  total  porosity  of  the  porous  medium 

0,  unsaturated  or  "air  filled"  porosity  of  the  medium 

X  hydrodynamic  dispersivity 

p  viscosity  of  permeate 

Pc  density  of  gas  phase  species  c 

Pg  total  density  of  gas  phase 

p,  density  of  the  solid  matrix 

Pl  density  of  the  liquid 

maximum  surface  contaminant  adsorption  level  (mass  per  unit  area) 
X  time  constant  for  complete  contaminant  adsorption 


Subscripts 
p  pore 

t  tree  whose  trunk  reaches  outer  edge  of  medium 

s  tree  whose  trunk  is  in  arbitrary  plane 

min  minimum  size 

max  maximum  size 
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